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Abstract 



We numerically analyze the transmission through a thin disordered wire 
of finite length attached to perfect leads, by making use of banded random 
Hamiltonian matrices. We compare the Landauer and the Thouless conduc- 
tances, and find that they are proportional to each other in the diffusive 
regime, while in the localized regime the Landauer conductance is approxi- 
mately proportional to the square of the Thouless one. Fluctuations of the 
Landauer conductance were also numerically computed; they are shown to 
slowly approach the theoretically predicted value. 
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Two main theoretical approaches to the study of electronic conduction in disordered wires 
at zero temperature have been developed, starting from ideas originally due to Thouless 
and Landauer, respectively. In the Thouless approach, the wire is regarded as a closed 
system, and its conductance is defined through the sensitivity of eigenvalues to changes of 
the boundary conditions In the Landauer approach |3|||] the wire is an open system, 

connected to perfect conductors, and its conductance is related to the transparency of the 
wire to electronic waves, via a formula, that can be derived from linear response theory 
||. In the metallic regime, the conductance is given by the Thouless formula, which can 
in turn be derived from the Kubo-Greenwood formula supplemented by certain random- 
matrix theoretical assumption ||. The relation between the above two approaches is a 
subtle theoretical problem. 

In this Letter we present a comparison of the Thouless and the Landauer conductance, 
on a model which allows for an efficient numerical computation of both. In this model, the 
wire is described by a 'Band Random Matrix' (BRM), that is, a real symmetric matrix H 
of rank L such that matrix elements 7^ only for \i — j\ < b+ 1, where b is the halfwidth 

2 

of the band. We chose Gaussian distribution for all non-zero elements, with variance ^- for 
those on the side diagonals and a 2 for those on the main diagonal . Strictly speaking, BRMs 
of this kind describe ID disordered wires with long-range hopping; however, they are also 
deemed to catch most of the general features of quasi-one dimensional models |7|. 

In the Landauer approach the wire is connected on both sides to ideal conductors, or 
"leads". In our model the electron dynamics in the leads is described by two semi-infinite 
BRMs with the same bandwidth b, with all matrix elements inside the band equal to 1. The 
Schrodinger equation is 

6+1 

^2 7~ti,i+k U i+k = EUi (1) 
fc=-(6+l) 

where Tiij are the non-zero elements of the Hamiltonian 7i, which is an infinite band matrix. 
Elements Tiij are the same as above defined hij for 1 < i, j < L and are equal to 1 otherwise. 
In this way, one can write 7i = H,q + V, where Tio is an infinite band matrix, with constant 

2 



unit elements inside the band, which describes free motion, and V is a scattering potential. 
The original BRM H is just the projection of the Hamiltonian "H on the subspace spanned 
by states 1 < i < L. 

Free waves (which are eigenfunctions of Tto) are of the form u m = -h=e tkm , with the wave 



/2tt 

number k obeying the dispersion law 

"MIT = (2) 

The left hand side of (^|) is a periodic function of /c of period 2n, which takes the value 
26+1 when its argument is a multiple of 2ir. For — 1 < E < +1 Eq.(|2|) has 2b solutions, 
which can be found analytically for E = 1, E = and E = — 1 or by numerical methods 
at any other value of E. For example, they are kj = |g- W ith j = 1, . . . , 2b for E = and 
fc s = 2m w ith s = 1, . . . , 6 - 1 and = ^±zl w i t h £ = o, . . . , b for E = 1. The various 
allowed values of the wave-number A; at a given energy E correspond to different scattering 
channels and in the range 1 > E > — 1 all scattering channels are open. We thus have a 
multi-channel scattering problem with exactly b ingoing and b outgoing channels, which, if 
1 << b << L , formally reproduces the case of a quasi-one dimensional wire. 

The scattering properties of the conductor are described by a unitary scattering matrix 
S, that relates incoming and outgoing amplitudes, Il{r) and Ol{r)'- 

Jr) \Or, 

where the subscripts L and R stand for the left and the right lead. The 5* matrix can be 
written as 

r t 
t' r' 

where t and t' are the transmission sub-matrices in the two opposite directions, and r and 
r' are the reflection sub-matrices. 
The Landauer conductance is M 

,2 
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£M 2 , (4) 



where tij are the elements of the t sub- matrix in (J3|). We have computed the S-matrix by 
standard methods, namely by numerically solving the Lippman-Schwinger equation 



u 



± - G± ■ Vu* = u (5) 



where u are free eigenfunctions, interacting ones. The free Green function, Gq — (E — 
Ho±ie)~l , has matrix elements 

{Gt)n, m = (n\(E - H ±iey 1 \m)\ e=0 



1 r-2-K 77 e i(m-n)k 
2^ h aK sin( 2i>+i fc) 

£±« — hS 7 



e=0 



which were numerically computed by exploiting the Residues Theorem. Finally, the scat- 
tering matrix was computed via the standard formula: 



Sij = 5ij - 2iri\ 



I dki 
dH 



I dkj 
dH 



<Ui\V\uf> (6) 

H=E 



H=E 

The unitarity condition S + S = I provides a check on the precision of the results, and was 
fulfilled within a typical error 1CT 10 . 

We have considered the dimensionless residual conductance 

2lTh r (7\ 
9 = —S (7) 



and investigated the scaling properties of the geometric average g av = exp (In g) in the regime 
1 ^ b ^ L, where L is the length of the disordered sample. The brackets (. . .) mean here 
average over different realizations of the potential (20000 realizations in our computations) 
while the bar means average over the energy range (—1,1). As expected, g av turned out 
to depend on b and L only through the localization parameter x — j- (Fig.l), where the 
physical meaning of b and b 2 is that they are proportional respectively to the mean free path 
and to the localization length. For x ^> 1, i.e. in the delocalized regime, this dependence has 
the form g av oc x oc L^ 1 , which is the "ohmic" behavior (Fig. 2). Instead, in the localized 
regime i C 1, g av ~ exp(— c/x) (Fig. 3). 

Our results can be summarized by the following interpolating law, valid in all regimes: 



a + frexp (— d/x) . , 

9av = j—j—, — (8) 

exp [c/x) — 1 

where c ~ 2.70 and d ~ 0.55 (obtained from the fitting in the localized region) and a ~ 6.58 
and b ~ —5.53 (obtained from the fitting in the metallic region). If the conductance is 
computed without performing the energy average, a similar scaling law is obtained, though 
with slightly different constants (curves a and b in Fig. 4). 

We now turn to the Thouless approach. In order to define the (dimensionless) Thouless 
conductance, the BRM which describes the disordered sample has to be 'periodicized' as 
described in ref. [f|. The 'periodic' BRM thus obtained describes a ring-shaped conductor, 
and depends on a phase <f> that has the meaning of an Aharonov-Bohm flux switched on 
through the ring. The disorder-averaged Thouless conductance is 



d 2 E 



(9) 

0=0/ 



where A is the mean level spacing, and (. . .) denotes statistical averaging. Although the 
scaling behavior of K, was described in ref. [B] by an interpolating law different from Eq. (^j), 
the latter law yields a slightly better fit with the parameters a' ~ 0.11, b' pa —3.1 • 10 -2 , 
d pa 1.56 and d! ~ 0.99. Therefore Eq. ([8]) can be used to describe the scaling of both the 
Landauer and the Thouless conductance, albeit with different parameters. 

A key problem in comparing the Thouless and the Landauer conductance is connected 
with the matching of the "free" and of the "disordered" part of the Hamiltonian, i.e., with 



the effect of "contacts", which was discussed in refs. |pTO| . If the arbitrary scale factor a 
appearing in the disordered part of our Hamiltonian TC is varied, while keeping the free 
Hamiltonian Tio unchanged, the Thouless conductance is not affected, nor is the value of 
the localization length in the 'infinite' sample. The Landauer conductance, instead, is to 
some extent modified. It becomes very small both at small and at large values of a, and 
it has a relatively broad maximum around a = 1. A value of a close to 1 is also found on 
requiring that the local density of states for the "free" dynamics be the same as the density 
of interacting states in the sample (having computed the latter from the interacting Green 



function, in a standard way). Our choice of a = 1, which maximizes the transmission, is 
very similar to the "matching wire" prescription, suggested by Economou and Soukoulis 
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The comparison between Landauer and Thouless conductance, computed as described 
above, is shown in Fig.4. In particular, in the metallic regime (x 3> 1) we have g av ~ (7.5 ± 
0A)K av +0.S, that is, the two conductances are proportional in the diffusive regime, where (||) 
is assumed to be valid. This difference cannot be removed by using of a different expression 
for the Landauer formula, as in ||. The meaning of the coefficent of proportionality ~ 



7.5 ± 0.4 is not clear to us ||12|| . It is of course possible that this coefficient becomes closer to 
1 in the limit b — > oo, L — > oo, b 2 /L = const.; this we were unable to check due to obvious 
numerical limitations, but, in the parameter range we were able to explore, the variation of 
this factor is quite slight. 

In the localized regime (x <C 1) we get \n(g av ) ps j3\n(K av ) + 5.66, that is to say, the 
Landauer conductance is proportional to K% v , with (3 = 1.7. The error in this numerical 
estimate is of order 0.1, due to the difficulty in computing curvatures in the localized regime. 
In this connection, we have also studied two different tridiagonal models: a BRM with 6 = 1 
and a standard Anderson model with Gaussian disorder on the main diagonal. Comparing 
Landauer conductances and curvatures on these two models we got Perm = 2.08 ±0.04 and 



ft And = 2.09 ± 0.01 fl3| . It seems therefore reasonable |L4| to conjecture (3 — 2. 

Let us now discuss the fluctuation properties of the conductance. First, we have found 
that in the localized region the following relation holds: Var (In g) = —K,(lng) + 7 with 
k = 2.00 ± 0.02, and 7 = —4.3 ± 0.5 (Fig.5). This has to be compared with the variance 
of the curvature, which in was found to obey a similar law, but with /c « 1. This fact 
reflects the most serious difference between Landauer conductance and curvature in the 
strongly localized regime. 

We have also performed a detailed analysis of the distribution P(g) of conductance in 
different regimes. Our results are in perfect agreement with theoretical expectations for 
quasi-lD conductors [£L5|], in both the localized and in the metallic regimes. In particular, 
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we have found that the shape of this distribution is controlled by the same scaling parameter 
x: it is log- normal in the localized region, x C 1, and normal in the delocalized one, x ^> 1. 

Special attention was paid to universal conductance fluctuations in the metallic regime. 
We have performed numerical simulations for large values of the scaling parameter ^, which 
is proportional to the ratio of the localization length (in the infinite sample) to the finite 
sample size. At the same time, we have taken the ratio ~ small, that is, we have taken sam- 
ples considerably larger than the mean free path. Our data for the variance of conductance 
versus the scaling parameter are summarized in Fig. 6. For not very large values of x ~ 5 
the value of Var(g) remains close to |, slightly different from the expected value ^. As 
discussed in |16| , the value | is predicted by RMT, while ^ follows from different theoretical 
approaches, based either on diagrammatic calculations, or on a diffusion-equation approach 
for transfer matrices; it is also obtained from RMT, if an appropriate repulsion law is used 
for eigenvalues of the transfer matrix. More detailed numerical studies, at larger values of 
x ~ 10, reveal that ^ is slowly approached on increasing x. All our results were carefully 
checked by using a very large number of realizations, and by implementing different random 
number generators. 

Our data indicate that the value ^ is approached very slowly, in such a way that actual 
data linger around | in a significant parameter range. That the approach to ^ can be 
relatively slow has been pointed out in Ref. ||10|| , where this effect was related to the coupling 
with channels. 



Finally we would like to remark that, since the velocity in the ith. channel is v j = ^ 



H=E' 



in our model different channels have different velocities. This raises the question, how does 
this fact affect our results. We have then modified the structure of the leads, in such a way 
that the dispersion law gives the same velocity in all channels. For this we have chosen the 
free hamiltonian with elements = dij + S i+ b + i t j + <5j_5_ij. The data obtained in this way 
are practically indistinguishable from those in Fig. 1. Therefore the fluctuation properties 
of conductance appear to be independent of the specific properties of the leads. 

In conclusion, we have studied the Landauer conductance for thin disorderd wires at- 



tached to perfect leads in ID (or in quasi-lD geometry). The model we have used is based 
on Band Random Matrices with the modifications imposed by the presence of the leads. The 
results are compared with those found for the Thouless conductance both in localized and 
metallic regimes. In the localized regime our results reveal that the Landauer conductance 
is proportional to the square of the Thouless conductance. In the metallic regime they are 
proportional to each other. Therefore, our data confirm theoretical expectations for the 
fluctuation properties of conductance in the metallic regime. 
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FIGURES 

FIG. 1. Logarithm of the geometric average of the Landauer conductance g av = exp(ln(g)) 
vs. logarithm of the scaling parameter x = ^. Numerical data correspond to different values of L 
and b in the ranges 50 < L < 1500, 7 < b < 80. Triangles are for the dispersion law (2), circles for 
a dispersion law giving the same velocity in all channels. The full line gives the fitting law eq.(8). 

FIG. 2. Geometric average of Landauer conductance vs. the scaling parameter ^ in the 
delocalized regime (b 2 » L). The linear fit shows the Ohmic behavior of Landauer conductance: 
g av = exp (ln(g)) = + B, with A = 0.44 ± 10" 2 , B = 0.6 ± 0.1. 

FIG. 3. ln(g av ) vs. the inverse of the scaling parameter in the localized regime (6 2 <C L). 
The linear fit shows the exponential decrease of Landauer conductance due to localization: 
(ln(#)} = -Cjp + D, with C = 2.70 ± 4 • 10~ 2 , D = 1.8 ± 0.2. 

FIG. 4. Comparison between the fittings of Landauer conductance, computed with (curve a) 
and without (curve b) an energy average in the interval (—1, 1). Notice that the effect of the energy 
average is very small. Curve c gives the fit of Thouless conductance computed in ref.[8]. 

FIG. 5. Variance of the logarithm of Landauer conductance g vs. < ln(g) >. The numerical 
fit shows that in the localized regime Var(ln(g)) = —K(ln(g)) + 7, with k = 2.00 ± 0.02 and 
7 = -4.3 ±0.5. 

FIG. 6. Variance of the logarithm of Landauer conductance vs. logarithm of the scaling 
parameter: the plateau in which Var(g) is independent of b corresponds to the regime of UCF, 
before the onset of the ballistic regime. The theoretical expected value ^ is reached only for 
x « 10, while for x s=s 5 data linger around | (upper and lower horizontal lines). 



11 




h-1 



' ^ — (Ni <N rn 



( AB §)ui 




<(S)up> 




o o o o o o o 

*sO in "st cn cn 



((S)ui)iba 



oo 



I u 



1 1 


i i 


1 


1 


M «• ■ 

i 


' M W-.-H. 

i 


■W---M- H-HH-H 

i 






i— H 


oo 




XT 


(N 






d 


o 


o 


o 


O 


d 


d 




d 


d 


d 


d 



O 

d 



